Steady-state many-body entanglement of hot reactive fermions 
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Entanglement is typically created via systematic intervention in the time evolution of an initially 
unentangled state, which can be achieved by coherent control, carefully tailored non-demolition 
measurements, or dissipation in the presence of properly engineered reservoirs. In this paper we 
show that two-component Fermi gases at ^ /iK temperatures naturally evolve, in the presence of 
reactive two-body collisions, into states with highly entangled (Dicke-type) spin wavefunctions. The 
entanglement is a steady-state property that emerges — without any intervention — from uncorrelated 
initial states, and could be used to improve the accuracy of spectroscopy in experiments with 
fermionic alkaline earth atoms or fermionic groundstate molecules. 



Many-body entangled states are known to be useful for 
quantum computing, quantum teleportation and cryp- 
tography protocols pQ , and precision metrology [2 . With 
these applications as motivation, the physics community 
has invested tremendous effort in preparing, stabilizing, 
and measuring entangled systems. Much of this effort has 
relied on coherent (Hamiltonian) dynamics to arrive at 
entangled states starting from less exotic states with only 
classical correlations. However, these approaches typi- 
cally suffer from the necessity to either carefully engineer 
interactions between particles or to prepare extremely 
pure and specific initial states (or both). A bottom up 
implementation of coherent control has yielded entangled 
states of up to 14 atoms with relatively high fidelity [3] (in 
ion experiments), and a top down approach has yielded 
weakly entangled states in a Bose Einstein condensate 
of ~ 10 4 neutral atoms [4]. A promising alternative to 
coherent control is the collective-nondemolition measure- 
ment of some observable with a finite variance in an ini- 
tially classical state. Such approaches have been used 
to generate entanglement (in the form of spin squeezing) 
amongst as many as 10 6 cold thermal atoms [5]. How- 
ever, collective and coherence preserving measurements 
are generically difficult to make, and the induced non- 
classical correlations are typically weak. 

In this paper, we show that two-component non- 
degenerate fermionic gases can be driven by reactive s- 
wave two-body collisions into steady-state spin configura- 
tions that, for a given value of the saturated particle num- 
ber, are pure and highly entangled. The entanglement 
comes in the form of Dicke states [6 , in which the spin- 
wavefunction is fully symmetric under interchange of the 
particles (with the burden of fermionic antisymmetry be- 
ing taken up entirely by the motional degrees of freedom). 
Such states have been sought in experiments with ultra- 
cold bosons for use in Heisenberg-limited phase measure- 
ments [7] , however these approaches typically suffer from 
the necessity to reach extremely cold temperatures (for 
the validity of a two-mode approximation in a double- 
well potential) or to employ Feshbach resonances [8J (to 
enhance spin exchange interactions for two-component 
Bose Einstein condensates). The only requirements to 



achieve such entanglement in the steady-state of lossy 
non- degenerate fermions are to have an SU(2) invariant 
single-particle Hamiltonian (in the pseudo-spin degrees 
of freedom) and a significant separation of timescales be- 
tween s- wave and p-wave collisions, with the second re- 
quirement typically being satisfied for temperatures in or 
below the fiK range. 

Because the desired property (i.e. the Dicke type of 
spin-entanglement) persists in the steady-state of dissi- 
pative dynamics, we do not rely on the highly controlled 
coherent manipulation that is typical of spin-squeezing 
experiments with bosons. Driven, dissipative prepara- 
tion of nontrivial steady-states has been considered be- 
fore in the context of many-body atomic systems [SHU] , 
and has been achieved recently in [12j [13] . In contrast to 
these examples, the mechanism described here is intrinsic 
and generic to a variety of interesting and experimentally 
relevant systems, such as fermionic alkaline-earth atoms 
(AEAs) and fermionic dipolar molecules, and does not 
require any special engineering of the system-reservoir 
coupling. After presenting calculations in support of our 
claims, we discuss the possible realization of such steady 
states in an experiment. In particular, we will propose a 
simple proof of principle experiment in which the steady- 
state entanglement can be revealed via Ramsey spec- 
troscopy of the 1 6'o to 3 Po clock transition of an AEA [14] . 
In this case, we will see that the interferometric precision 
stays relatively constant even as most of the particles are 
lost (all but ~ ^/J\f in the long time limit), signaling the 
development of quantum correlations and the pursuant 
violation of the standard quantum-limit. The total loss 
of precision (due to loss of particles) exactly cancels the 
gain due to the growth of quantum correlations. How- 
ever, a persistent precision under loss of particles can 
provide enhanced spectroscopic accuracy; in particular, 
the steady-state under consideration is largely devoid of 
mean- field clock shifts. 

Our description of spin-^ fermions with two-body re- 
active collisions relies on the formalism detailed in Refs. 
[T5UT7] . generalized for fermions, where we assume the 
temperature to be sufficiently low that losses are domi- 
nantly in the s-wave channel. As in Ref. [15 , large ki- 
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netic energy of fermions in the outgoing channels (which 
for reactive molecules can correspond to temperatures 
in the 10K range) guarantees they will be rapidly lost 
from any typical atom trap, justifying a Born Markov 
approximation. Given a density matrix g for the sys- 
tem (fermions, Hilbert space 5?) plus reservoir (outgoing 
channels of the inelastic collisions, Hilbert space ^), the 
Born-Markov approximation leads to a master equation 
for the system reduced density matrix p = Ti^[g] [16] : 
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ApM~\ Jd 3 r(j^Jp + pJ^J-2JpJ^). (1) 

The system Hamiltonian H = Ho + g J d?rj^ J is com- 
posed of an unspecified single-particle term Ho and an 
interaction term with coupling constant g = Anh 2 aR/m, 
where m is the particle mass and a = clr + iaj (aj < 0) 
is the complex s-wave scattering length. The jump oper- 
ators are defined by J(r) = i/j^(r)i/j^(r) (their explicit r 
dependence is suppressed in the integrals above), where 
^cr(r) annihilates a fermion located at position r in in- 
ternal state a G {t?l}> an d ft = —^h 2 ai/m. Assuming 
without loss of generality that the initial number of par- 
ticles J\f is even, the relevant system Hilbert space can 
be written as a direct sum over spaces with well defined 
particle number, = y M © y M ~ 2 © • • • © S*°, between 
which coherence never develops. Hence, the density ma- 
trix can be decomposed into a sum of density matrices 
in each particle-number sector, any one of which we la- 
bel by p n once normalized. Furthermore, any Hilbert 
space Sf n can be decomposed into a direct product be- 
tween motional (m) and spin (s) degrees of freedom, 
cfn _ ^ ^ anc j we can define a reduced spin den- 
sity matrix by p™ = Tr y m [p n }. For what follows, it will 
be useful to define a fidelity in a given Dicke state [28] of 
the spin degrees of freedom of n particles, \S = n/2, S z ), 
given by the population of p™ in the Dicke state 
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Here S and S z are quantum numbers for the total spin 
and its projection along the z-axis, respectively. 

Two particles. To make the physics clear in a simple 
context, we begin by considering two fermions in a single 
double well potential (which could be formed in an opti- 
cal super-lattice p~8j[19]). We consider a single wavefunc- 
tion ip a (r) in each well (a G {L, i?}), denote the creation 
operator for a fermion in spin state a and wavefunction 
(p a by V4qo and choose an initial state i/j^ L ^ R \vac) with- 
out spin correlations. Within a tight binding model for 
these two wavefunctions, the Hamiltonian is 

cr a=L,R 

and the master equation reads 

a=L,R 

(4) 
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FIG. 1: (Color online). The solid red line is the fidelity of 
the final density matrix (after post-selection for a non- vacant 
well) with respect to the S z = Dicke state, and approaches 
1 (black dotted line) at times long compared to 7 -1 . The 
oscillations imposed over the exponential decay are due to 
the inter- well hopping. 



Here J is the inter- well hopping, J a = ^^ a ^^ a1 U = 
g j d?r\ip a (r)\ A is the onsite interaction energy, and 7 = 
n J d 3 r|(^ a (r)| 4 is the onsite loss rate. The initial state 
can be decomposed into an evenly weighted superposition 
of triplet and singlet [(^l L ^t^±^{ L ^l^)|mc), with plus 
for the triplet], and the spin wave function of the triplet 
is the entangled Dicke state |1, 0). The triplet, having 
a spin wavefunction that is symmetric under exchange, 
has an orbital wavefunction that is antisymmetric under 
exchange, and therefore it is "dark" to s- wave losses, by 
which we mean simply that Jl\^) = Jn\t) = 0. It also 
happens to be an eigenstate of and so it is stationary 
under propagation of the master equation Q. On the 
other hand, there are no dark eigenstates in the singlet 
sector, and as a result p 2 s is pure at long times and sat- 
isfies Ti,o = 1. In other words the steady state of our 
system, when restricted to the subspace with two par- 
ticles, is the entangled Dicke state |1,0). It should be 
noted that there is also a 50% probability of obtaining 
the vacuum, and hence in an array of double wells the 
entanglement fidelity is only unity after post-selection 
of the non- vacant wells. In this simple example we see 
an important general feature of the physics we will dis- 
cuss, that even purely local (intra-well) dissipation, when 
coexisting with Hamiltonian dynamics that delocalizes 
the particles, generates non-local (inter-well) spin corre- 
lations in the steady state. 

Many particles. Solving Eq. ([!]) exactly for initial 
states with AT > 2 initial particles quickly becomes 
impossible, but strong statements can nevertheless be 
made regarding the steady state. It is crucial to ob- 
serve that the jump operators only remove spin singlets 
from the system, which follows from Fermi statistics com- 
bined with the even exchange symmetry of the spatial 
part of any two-particle wave function susceptible to s- 
wave scattering. Intuitively, this suggests that losses 
do not decrease the expectation value of the total spin, 
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& = \ I d 3 rip^(r)T aa 'ip a ,(r) ( T being a vector whose 
components are the Pauli matrices) . Mathematically, we 
say that ^(S • S) = Tr[pS • S] = 0, which can easily be 
verified in the case when T~L is SU(2) invariant by check- 
ing that [S • £, J{r)] = 0. A stronger consequence of the 
commutation of all J{r) with S • S is that population 
in any sector of total spin, Vs, is also conserved. Be- 
cause any state with well defined total spin S must have 
(TV) > 2S particles (where N = J d 3 r^ip . is the to- 
tal number operator), an immediate consequence is that 
the loss of particles can only yield the vacuum deter- 
ministically at long times if the initial state is a total 
spin-singlet. For an uncorrelated spin state, such as a 
non-degenerate thermal distribution of J\f fermions in a 
balanced incoherent mixture of f and ^, it can be shown 
that [20] 

^ Tr^r \^ + ll 

N(t) ee Tr[pN] > MPs = ^fttr ] ~ ^ 

which places a lower bound on the steady-state expecta- 
tion value for the number of particles N(t). This expec- 
tation value determines the particle number in a typical 
steady-state configuration, and is achieved (on average) 
without any post selection, but variations of the steady- 
state particle number will occur from shot to shot. Tak- 
ing Stirling's approximation for large M yields an approx- 
imate bound N(t) > ^ttJV/2. For the chosen restriction 
on the initial state, the validity of Eq. ([5| depends only 
on the SU(2) invariance of and not on its precise form. 
Whether the bound ([5| is saturated in the steady-state, 
however, is an important and delicate issue; an affirma- 
tive answer guaranties that all of the describe pure 
Dicke states in the steady-state. Demonstrating that this 
bound is indeed saturated in certain experimentally rel- 
evant situations, namely a ID harmonic trap and a ID 
Hubbard chain (optical lattice), is a central technical re- 
sult of this paper. 

Saturation of the bound in Eq. ([5| is guaranteed if, 
for any fixed value of n and S Z1 the pure density matrix 
|n/2, S z )(n/2, S z \ is the unique steady-state reduced spin 
density matrix. This uniqueness, in turn, is equivalent to 
requiring that any dark-state with quantum numbers n 
and S z has a well defined spin wavefunction given by the 
Dicke state \n/2,S z ). In the supplement we prove this 
to be true for a ID harmonic oscillator potential, and 
we have verified it numerically for a ID Hubbard chain 
(see below). It is worth noting at this point that, while 
the equivalence of dark-states with the Dicke states is 
intuitive, there are natural Hamiltonians for which this 
intuition is incorrect. In particular, all Hamiltonians in 
D > 1 that are separable in cartesian coordinates do have 
dark-states with J r s,s z < 1- 

In order to verify the above statements numerically, 
we have performed quantum trajectories simulations for 
an 8-site Hubbard chain with open boundary conditions, 
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FIG. 2: (Color online). Calculation of particle number [N(t), 
solid red line] and average Dicke state fidelity [T{t), dashed 
blue line] for an 8-site Hubbard chain via quantum trajec- 
tories. For the number, the shaded region is an estimate of 
the statistical error from sampling of a finite number of tra- 
jectories. The black dotted line is the analytic bound in Eq. 

an initial filling of one particle per site, and zero po- 
larization (AT = 8 and S z = 0). In general we have 
lower-densities in mind for any experimental application, 
but using one particle per site allows us to stretch the 
numerics to the largest M possible. In Fig. [2] we show 
the calculated particle number and average Dicke state 
fidelity, F{t) = \Y^s=i^s,o^ an d one can see that the 
former saturates the bound Eq. ([5| while the latter ap- 
proach unity at long times. For this calculation we solve 
for O(10 4 ) trajectories with no approximation. 

Experimental realization. Dicke states are known to 
be useful for a variety of quantum information proto- 
cols, including but not limited to quantum secret shar- 
ing [21], teleportation [22], and sub shot-noise limited 
precision spectroscopy [7]. Here we give a brief descrip- 
tion of how the proposed Dicke state preparation could 
be used in precision spectroscopy of the clock transition 
in alkaline-earth atoms. For a fixed interrogation time, 
spectroscopy on AT uncorrelated atoms has a phase sen- 
sitivity Sep > I/VAT, a bound known as the standard 
quantum limit (SQL). This bound can be understood as 
the minimum tipping angle needed to cause a coherent 
spin state (CSS) to have an uncertainty cone that pre- 
cludes its initial position [23] (to one standard deviation). 
On the other hand, spectroscopy on a Dicke state of M 
particles with spin S z = has the potential to reach the 
Heisenberg limit (HL) of phase sensitivity, Sep ~ 1/Af 
[3 [23]. It is important to realize that the production of 
Dicke states with \fj\f fermions via two-particle loss does 
not actually enhance the phase sensitivity relative to the 
initial state with AT fermions; the enhancement in phase 
sensitivity between the SQL and HL exactly compensates 
the reduction of particle number. However, the reduced 
particle number in the Dicke state and darkness to real s- 
wave interactions (which if present generate clock shifts) , 
can render the accuracy of the final Dicke state superior 
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to that of the initial J\f fermion uncorrelated state. 

Rather than allowing losses amongst a macroscopic 
sample of atoms, for which the approach to the steady 
state could be quite slow, we imagine an array of T ID 
tubes created by a 2D optical lattice. Although there will 
be variations in the atom number from tube to tube, for 
simplicity we take each tube to have exactly J\f fermionic 
AEAs in the 1 So electronic state and I z = I nuclear-spin 
state, denoted l 1 ^,/). For the analysis in this paper to 
be valid, the temperature should be small compared to 
the vibrational level-spacing in the transverse tube direc- 
tion, and also low enough that only the harmonic part of 
the trapping potential along the tube axis is sampled by 
the atoms. A 7r/2-pulse on the spin degrees of freedom 
[I 1 S , 1) -^(|%,/) + |%,/-l>)], followed by single 
particle dephasing [29] , generates a statistical mixture of 
the two spin states (I z = I and I z = I — 1). Losses 
can be initiated by applying a 7r-pulse on the clock tran- 
sition (^So,/*) \ 3 PoJ z }). We estimate that this tt- 
pulse [30 can be achieved on the < 100/is timescale with- 
out exciting transverse excitations in the tubes (which, 
if present, violate the assumption of a ID geometry and 
destroy the uniqueness of the steady-state). Thus the 
transfer into 3 Po is sufficiently fast that it can be con- 
sidered instantaneous on the initial timescale of reactive 
collisions — which, based on universal considerations for a 
Lieb-Liniger gas, we estimate to be > 1ms for experimen- 
tally relevant ID densities [15] — such that it suddenly 
initiates strong 2-body s-wave losses. 

The steady state of the system is a statistical mixture 
of Dicke states in the different tubes, each having some 
value of Vj particles (centered around Vq & \fj\f) and 
spin projection Sj (centered around zero). Spin selective 
transfer of \ 3 Pq,I — 1) into I 1 So, I) maps the spin degree 
of freedom onto the clock states, leaving a spin-polarized 
sample, and Ramsey spectroscopy on the clock transition 
can then be performed [23]. Despite the fluctuation of 
both Vj and 5| from one tube to another, it can be 
shown (see the supplement) that the minimum resolvable 
rotation angle in a Ramsey experiment scales as 



1/VoVf. 



(6) 



This result can be interpreted as the existence of Heisen- 
berg limited sensitivity for each tube, which is then com- 
bined between tubes in a statistically independent man- 
ner (hence the 1/y/T). In order to utilize this phase 
sensitivity, the initial value of S z = S z for the entire 
ensemble must be accurately known. Because S z is con- 
served by the losses, it can be measured before transfer 
to the 3 Po, and hence the measurement does not need to 
preserve any inter-particle correlations (since these de- 
velop during the losses). Accurate measurements of this 
type and precision for ~ 100 atoms in an optical cavity 
have recently been demonstrated [24] . 

The primary limitations on the final state fidelity 
achievable in experiments is likely to be a combination 




FIG. 3: (Color online), (a) An array of T ID tubes, each 
having Vj atoms in a Dicke state, (b) Bloch sphere represen- 
tation of a Dicke state in a particular tube. 



of finite p-wave losses (which the Dicke states are not 
dark to) and magnetic field gradients. At sub /iK tem- 
peratures, the s-wave losses in a spin mixture of 87 Sr are 
expected to be about an order of magnitude faster than 
the p- wave losses [25 . For reactive molecules (or 171 Yb), 
where the inelastic collisions are expected to more fully 
saturate the unitarity bound [26] , this separation of rates 
will most likely be even larger. Magnetic field gradients 
couple sectors of different total 5, all of which are sepa- 
rated from the Dicke manifold by a gap for finite systems 
and nonzero or, so in principle their adverse effects can 
be suppressed to first order [27]. Furthermore, if the two 
components of the Fermi gas are two nuclear spin states 
of an AEA, they will be extremely insensitive to mag- 
netic field gradients: We estimate that typical gradients 
(lmG/cm) will cause spin dephasing on a 100s timescale 
for a linear system size of 100/im. This timescale is sev- 
eral orders of magnitude longer than the initial two-body 
loss rate in tightly confined ID tubes, which we estimate 
to be on the order of 10ms for 87 Sr (assuming a 50£?r 
2D lattice and scaling the density dependent loss rate 
from Ref. [25]), and even faster for 171 Yb [26]. A more 
quantitative analysis of the effects of both magnetic field 
imperfections and finite p-wave losses requires numerical 
simulations beyond the scope of this work, and is left for 
future study. 

Conclusions. In this paper we have demonstrated 
that fairly unrestrictive initial conditions, without in- 
tervention or engineering and in the presence of re- 
active two-body collisions, are sufficient to generate 
steady-state spin entanglement between non-degenerate 
fermions. These reactive collisions, which occur both in 
optically excited alkaline earth atoms and many dipolar 
molecules (e.g. KRb), are typically viewed as an imped- 
iment to interesting physics, but clearly this need not be 
the case. We expect this physics to enable the distilla- 
tion of Dicke states from initially uncorrelated fermionic 
atoms and molecules, hence extending the scope of a va- 
riety of experimental progress made in the spin squeezing 
of bosons. 
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Uniqueness of the steady state 

Here we show that for a ID harmonic oscillator, in a 
particular sector of Hilbert space and for a particular 
value of S z , the unique steady-state reduced spin density 
matrix is given by 

p n 8 =\n/2,S')(n/2,S'\. (SI) 

The extension of what follows to the ID Hubbard chain is 
fairly straightforward, and will be described in more de- 
tail in future work. As discussed in the text, it is sufficient 
to prove that all dark eigenstates of the non-interacting 
Hamiltonian 

Hho = Jdx ^t(x) + l -muj 2 x 2 ^j ^(x) (S2) 

have a maximally symmetric spin wavefunction. To un- 
derstand the properties of its dark eigenstates under par- 
ticle exchange, we will actually work in first quantization 
writing an eigenstate for M particles as 

* = Y, A 3**(ri,---,w)\S). (S3) 

a 

Here the j th component of the vector <j, aj G {t?l}? 
labels the spin orientation of the j th particle (along some 
arbitrary quantization axis, which we'll call z), and the 
total spin wavefunction in any term of the sum is 

|<?) = |cri> <g) |cr 2 > <8) - - - (8) |o>/>. (S4) 

The sum over a should be understood as independent 
summations over each index 

E = EE-E- (ss) 

a <J\ CT2 (7j\f 
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the coefficients Aa are arbitrary, and <I>^ is a normalized 
orbital wavefunction for the M particles. Dark states of s- 
wave losses have zero expectation value in the interaction 
operator 

u = Urnn = g Yl s ( rm rn )' ( S6 ) 

m<n m<n 

and this expectation value can be evaluated as 
U = J @r^*U^ 

= El^-I 2 / yr^UGt, (S7) 

a 

with Q)r = Ylj drj . The last equality holds because the 
interaction is spin-independent. It is crucial to realize 
that the operator U is posit ive-semidefinite, which means 
that satisfying U = actually implies the stricter con- 
straint 

j @r<$>* a U$z = W. (S8) 

In addition to the operator U being positive-semidefinite, 
the constituent pairwise interaction operators are as well. 
Hence, the condition J @r QtUQj = actually implies 
that 

J ®r <S>%U mn $t = Vm/ n. (S9) 

The above set of equalities can now be used to pin down 
properties regarding the exchange symmetry of the wave 
functions <I>^. 

Implications for the exchange symmetry of the orbital wave 
functions 

Let's choose two particles, say particle 1 and particle 
2, and define relative and center-of-mass coordinates for 
them as r = ri — r2 and R = (r\ + r 2 )/2. We can then 
expand the wavefunction <1>^ in a way that makes the 
symmetrization with respect to exchange of particles 1 
and 2 explicit: 

§ B = ^2Bt(3 S x (p a (p i3 ) (r)(pp(R)ip s (r 3 , ...,r N ). (S10) 

(3s 

In the above (p a are harmonic oscillator wavefunctions, 
the ip s are a complete and orthonormal set of eigenfunc- 
tions for the remaining M — 2 particles, and the notation 
a(/3, s) implies that the relative wavefunction of particles 
1 and 2 is uniquely determined by j3 and s. This point is 
crucial, and relies on the observation that for ^ to be an 
eigenstate, when expanded in terms of eigenstates & a all 
of the eigenstates must have the same eigenvalue. Hence 
the energies of states a, /?, and s are constrained to add 



to some fixed value. From now on we'll drop this explicit 
dependence. The interaction energy between particles 1 
and 2 is given by 

U 12 = J @r& a U 12 *a 

= ^fi ff> ^ Q (0)xfi^X(0) 

= J2\B^ a (0)\\ (Sll) 

(3s 

and hence Ui 2 = implies that 

B<?p a <P«(0) = V/3,5 (S12) 

So B#p s must be zero for all even wave functions (all of 
which are finite at the origin) , implying that is strictly 
odd under interchange of particles 1 and 2. By repeating 
the above argument for two arbitrary particles m and n, 
it is easy to see that is strictly odd under interchange 
of any two particles. 

Implications for the spin wavefunction 

We can now ask what the antisymmetry of implies 
for the full wave function 

9 = ^Aff*ff(ru...,r N )\ff). (S13) 

Under interchange of two arbitrary particles we have 
\a) — » | a'), and we obtain the new wave function 

*' = -Y.As$>s{ru---,r N )\a') 

a 

= ~Yj A ^ ®ff'(ri,...,r N )\a) 

a 

= -Y J A s $ s (r u ...,r N )\a), (S14) 

(7 

implying that 

As ®ff(r u ...,r N ) =Ab> ®ff>(r u . . . ,r N ). (S15) 

The second equality follows because switching a <-> a' 
in the summand just changes the order of the terms in 
the sum, and the third equality follows from the anti- 
symmetry of the total wavefunction under particle ex- 
change. By repeated permutations of various particles, 
this chain of logic can be used to demonstrate that all of 
the Ab &a( r i, • • • ? t n) are equal, and hence we have 

9=A*(r 1 ,...,r N )^\&). (S16) 

a 

Now A is just some normalization, which is related to the 
total z projection of the spin, and it is clear that ^ breaks 
up into the product of a completely antisymmetric orbital 
wavefunction multiplied by a completely symmetric spin 
wavefunction. 
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Phase sensitivity of the steady-state 

In order to estimate the phase measurement sensitivity 
of an array of ID tubes in the steady state, we begin by 
considering just the j th tube, with initial particle number 
A/", final particle number T>j after relaxing to steady state 
via collisional loss, and initial (and final) spin projection 
Sj. The initial M atoms can most easily be prepared in 
an incoherent mixture of spin up and spin down by sim- 
ply allowing a coherent state initially prepared along the 
x-direction to undergo single particle dephasing (which 
could be briefly enhanced via a myriad of methods). For 
the coherent state, the probability of a given 5J is given 
by a binomial distribution, which for large M is approx- 
imated by the continuous probability distribution 

V (S*) = yX e -2(sj/v^ ( si7) 

Note that without dephasing into a mixture, a coherent 
state of fermions does not undergo s-wave collisions. In 
the sense that such a distribution is easily prepared ex- 
perimentally, we take this to be a worst-case scenario; a 
distribution of Sj more sharply peaked around Sj = 
will enhance the phase sensitivity. The probability distri- 
bution of steady-state particle numbers in the j th tube, 
conditioned on a particular value of SJ, is given by 




(S18) 

For Sj = 0, this distribution is peaked around Vq ~ 
(giving the expected value of n(oo) quoted in the text). 



For \Sj \ > 0, the step function 6 reflects the fact that 
as particles are lost (remember that SJ is conserved by 
the losses), the remaining particles are maximally spin 
polarized once Vj = 2|SJ|. The second exponential pro- 
vides the proper normalization, \ J °° dVj V(Vj\S z ) = 1, 
where the factor of | comes from converting sums into 
integrals while respecting our assumption of even particle 
number. 

For this single tube, small rotations about the x-axis 
by an angle Sep cause a standard deviation in the final 
distribution of Sj given by [23] 

v(5<p, Vj , SJ) « 5<p^ ^ ^ ^ . (S19) 

for large Vj. For Sj = (before the rotation), this 
demonstrates that a discrepancy in Sj (after the rota- 
tion) of order unity is expected for Sep ~ hence the 
Heisenberg limited phase-sensitivity within a single tube. 

Our estimation of the phase sensitivity for an array 
of tubes relies only on the assumption that the initial 
(i.e. before the losses) value of S z = ■ Sj is known 
to within an uncertainty E, but does not require any 
knowledge of Sj in the individual tubes, greatly relaxing 
the experimental requirements. This uncertainty E guar- 
antees that, in principle, rotations causing deviations in 
S z of order E can be detected. For T tubes with well 
defined (i.e. measured) total spin projection S z , the ex- 
pected standard deviation in S z (denoted cr tot ) due to a 
rotation by angle tp about the x-axis satisfies 



T times 



T times 



(S20) 



Here 6 is the Dirac <5-function, reflecting the correlations 
established between the various Sj by the knowledge of 
S z , and the factors of | again come from converting sums 
into integrals while respecting the assumption of even 
particle number. If we ignore this i5-function constraint, 
which is valid in the large T limit, the integral simplifies 
greatly to 



r 



/oo roc 
dS* / dVjVCDjlS^ViSj) a(5tp, T>j,Sj) 2 
-oo ^0 



(S21) 



Equation (S21) can be evaluated explicitly to reveal 



4<x t 2 ot « TAT5(p 2 « TV 2 5v 2 



(S22) 



Setting the total standard deviation to E gives a mini- 
mum phase sensitivity of 



5cp n 



2E 



(S23) 



As has been demonstrated recently in Ref. [24], E ~ 1 
is possible for ~ 100 atoms in an optical cavity, as long 
as measurements that do not preserve coherence between 
the atoms are acceptable. Such measurements certainly 
are acceptable before the losses take place, since we re- 
quire no inter-particle correlations in the initial state 
(they develop dynamically due to the losses). There- 
fore, as quoted in the manuscript, we expect a minimum 



phase sensitivity of S(p m [ n ~ ~ ^/j" ^ achievable in experiment. 



